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We show that the density profile of a Fermi gas in rapidly 
rotating potential will develop prominent features reflecting 
the underlying Landau level like energy spectrum. Depend- 
ing on the aspect ratio of the trap, these features can be a 
sequence of ellipsoidal volumes or a sequence of quantized 
steps. 

Currently, there is an intense effort to cool trapped 
Fermi gases down to the degenerate limit. Recent ex- 
periments at JILA on 40 K has reached one half of its 
Fermi temperature 0. One of the motivations of cool- 
ing Fermions is to reveal their possible superfluid ground 
states, which can be quite novel in the case of multi- 
component Fermion systems Qj. Current theoretical es- 
timates, however, indicate that the interactions between 
different spin states of 40 K are all positive Q, implying a 
normal instead of superfluid ground state. The absence 
of superfluid ground states, however, does not mean that 
the system can not have novel macroscopic quantum phe- 
nomena. Quantum Hall effect is an excellent example. 
In a strong magnetic field, the energy levels of a two 
dimensional electron system organize into highly degen- 
erate Landau levels, leading to a whole host of dramatic 
effects. In this paper, we show that similar organiza- 
tions will take place in a fast rotating three dimensional 
trapped Fermi gas, leading to many macroscopic quan- 
tum phenomena. 

For neutral atoms in rotating harmonic traps, we shall 
see that Landau level like energy spectrum will appear 
when the rotation frequency f2 approaches the transverse 
confining frequency luj_ . This "fast rotating" limit might 
appear hard to achieve as the system is at the verge of fly- 
ing apart due to centrifugal instability B . Such instabil- 
ity, however, can be prevented by imposing an additional 
repulsive potential which dominates over the centrifugal 
force beyond certain radius. With centrifugal instability 
eliminated, the Landau-like levels will show up in many 
ways. We shall see that for cylindrical harmonic traps 
with u) z « luj_, where luj_ and uj z are the transverse and 
longitudinal trapping frequencies, the density is a sum 
of one-dimensional like density distributions residing in 
different "Landau volumes" . For traps with lo z compa- 
rable or smaller than u/j_, the density consists of a set 
discs along z, each of which is made up of a sequence of 
density steps quantized in units of Muj± / (jrh) . 

2D case: We first consider 2D rotating Fermi gases 
in harmonic potentials since they illustrate the basic 
physics. The Hamiltonian in the rotating frame is 



h-ql z = -^pi 



^Mu)\r 2 - fiz ■ r x p ± , (1) 



and eigenvalues of eq.(0) are 

e \w\ 2 /2gmQn e ~\w\ 2 

U n ,m(r,9) = 



(2) 



(3) 



£n,m = h(u± + fl)n + h(u>± — 0)m + LU±, 

where n, m are non-negative integers 0, 1,2,..; w = (x + 
iy)/a±; a± = y / h/Moj± ; and d± = (a±/2)(d x ±id y ). To 
derive eqs.(^) and (|^), we note that eq.(jl]) can be written 
as n 2 /2Af+(wj_ — U)L Z with n = pj_ — Mu>j_zxr, which 
is precisely the canonical momentum II = pj^ — x r of 
an electron in a magnetic field B in the symmetric gauge, 
with eB/Mc = 2ui±. The eigenfunctions of II 2 /2M are 
those in eq.(||) ||, with eigenvalues e^ m = huj±(2n + 1), 
where n is the Landau level index, and m is an "angular 
momentum" index labelling the degeneracy in each level. 
Since L z u n , m — h(m—n)u n ,m, eq.(|2j) is also an eigenstate 
of eq.(0) with eigenvalues eq.(|3|). Note that the function 
Un,m in c q-(0) peaks at 



= a 2 L {7i + m + 1), 



(4) 



where = (p x ,Py) and r = (x,y). The eigenfunctions 



and decays away as a Gaussian over a distance aj_. 

Eq. ([|) shows that the system is unbounded when Q > 
uj± unless an additional repulsive potential K, a ii(^) (say, 
introduced by an additional optical trap) is present. We 
shall in particular consider potentials Vwaii(r) which are 
zero for r < R but become strongly repulsive for r > R, 
with R >> a±. The specific form of Vw a n is not impor- 
tant for the key features discussed below, as long as it is 
smooth over length scale a±. The condition R » aj_, 
however, allows us to fit many m states inside r < R and 
is a necessary feature for many effects discussed below. 
Since V wa \i(r) is cylindrically symmetric, the eigenstates 
are still labeled by quantum numbers (n,m). For states 
originally with r n>m < R, eqs.(^|) and (g) remain valid 
because V wa \\ = for r < R. For states (n, m) origi- 
nally peaked beyond R, their energies increase rapidly be- 
cause V wa \\ is strongly repulsive. (For 40 K in a tight trap 
lu± = 4000Hz and 10 5 Hz, we have a± w 2.5 x 10~ 5 cm 
and 5 x 10 _6 cm resp. The condition R >> a± is satisfied 
for R > 5 x 10- 4 cm.) 

Let us first consider the case Q < lo±_ with Vvaii = 
0. The density in the ground state is p(r) = 
E„, m \u n ,7n(r)\ 2 &(fi-e n , m ), where 8(x) = 1 or if x > 
or < 0, [i is the chemical potential related to the parti- 
cle number N as N = J2 n m 0(/i — e„, m ). We can write 

P( r ) = J2n=o Pn( r > m *n)' where p n {r\L) is density contri- 
bution of the n-th Landau level with angular momentum 
states filled up to in = L: 



1 



Pn(r;L) 



L 

E 

m=0 



\u n , m (r)\ 2 Q(p - e n m ), 



(5) 



to* is the highest angular momentum state in the n-th 
Landau level with energy less than p, and n* is the high- 
est Landau level below p, 



to* = Int 



p/Ti — ujj_ — (lu± + Cl)n 



(LUj_ - CI) 



(6) 



Int 



where Int[x] denotes the integer part 

of x, and x is understood to be positive. Since (n, to* ) is 
the state in p n farthest from the origin, its peak location 
(r n = r n,m* n = a±\M + to* + 1) gives the size of p n . 
When Cl is very close to uj±, we have to* >> 1 and 



p- hCl(2n + 1) 
Mwi(wi - fi) 



(7) 



Note that the difference in area between successive Lan- 
dau discs is a constant 



7r(r 



Using eq. 



rl) = (Troi) 



2S1 



it is straightforward to show that 
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2na 2 ± 



1 — erf 



(1 



(8) 



(9) 



where s = (r/aj_) 2 — TOq, erf(x) = (2/y/W) J* e~ z dz, and 
the term [..] in cq.(^J) is of order (12toq) _1 and smaller. 
The densities p n (r) of the higher Landau levels (n > 0) 
can be generated from p as Q / o„(r;m*)= ^ (jV^)" 

[po^— a ±A; m* ) e A2 ]^ =Q . ,/From the properties of erf (x), 
it is clear that po is a constant l/(7ra^) within a disc 
of radius ro = ^/to.q + laj_ and has an edge of width 
A ~ "7f a -L [§)• If "^o >> 1' ^0 << r , and p can be 
approximated as a step function on scales larger than Ao. 
Likewise, p n can be approximated as a step function with 
somewhat larger width within the same approximation. 
Thus, when to* >> 1, we have 



(10) 



Eq.(|10|) and (|S|) then imply p(r) = J ^I{r), with 
7(r) = 9 (p - Mwi(wi - fi)r 2 - 7ifi(2n + 1)). Us- 
ing the identity Int[ir + 1] = J2 n ®{ot{x—n)), for all x > 
and a > 0, p(r) can be simplified to 



p{r) 



Int 



Ai - Mwifwi - f2)r 2 + ?iSl 



2M] 



(11) 



It is, however, instructive to re-derive eq.(pT|) in a way 
generalizable to arbitrary potentials. We rewrite eq.(Q) 



H — ClL, 



(pi — MClz x r) 2 1, 9 9s 9 . , 
2M L + 2 M ^- n >" (12) 



The first term gives a set of Landau levels with spacing 
2hCl, each of which contributes (ttA 2 )^ 1 to the density, 
where A 2 = h/ [MCI). If the second term in eq.(|lj) were 
absent, the density is given by p = I/(nA 2 ), where / 
is the total number of Landau levels below the chemical 
potential, / = Intf^jp]. When Cl is close to the 
second term in eq.(12) is slowly varying over the scale of 
A ~ aj_ and can be absorbed in the chemical potential. 
The density profile within local density approximation 
(LDA) is then 



p(r) 



MCI 
itTi 



J(r), I(r) = Int 



p(r) + hCl 



2HCI 



(13) 



where p(r) = p — ^M(uij_ — Cl 2 )r 2 . Clearly, eq.jl^) is 
equivalent to (O) up to correction (1 — Cl/u>±)<< 1 as 
Cl — ► ujj_. Eq.Ql3|), once established, is easily generalized 
to other potentials. In the presence of V wa n, one simply 
replaces p{r) in eq.([TJ|) by 

p{r) = p- ^M(usl - Cl 2 )r 2 - U wal i(r). (2D) (14) 

Eq.(|l3|) and eq.(|IJ) constitute the LDA solution for the 
2D rotating Fermi gas for both Cl < u>± and Cl > cu±. 
The schematics of LDA is shown in fig. 1(a) and 1(b). 

To understand the validity of LDA (eq. (ft3|V) , we have 
calculated the density numerically using eqT(g). The re- 
sult for a system of 2000 Fermions at Cl/ui± = 0.996 
is shown in figure 2a. The system exhibits a sequence 
of quantized steps at locations well described by LDA. 
The evolution of the density within the range 0.99 < 
Q/gj± < 1 of this Fermion system is shown in fig. 2b. As 
Cl decreases, more Landau levels are populated while the 
steps near the surfaces are closer together, (as expected 
from eq.(^) and (^|)), yet the step structures remain dis- 
cernable and correctly described by LDA, (the LDA con- 
struction is not displayed so as to keep the fig. 2b read- 
able). The behaviors of the densities at lower frequency 
0.98 < CI/ujx < 1 are shown in fig. 3a and 3b for a sys- 
tem of 1000 Fermions. At the lowest frequency displayed, 
(fig. 3a), the step structure near the surface is completely 
smeared out by the spread of the edges. Nevertheless, 
the density of the innermost plateau remained quantized, 
with a size correctly described by LDA. This core of quan- 
tized density (or "quantized core" for short) is a clear evi- 
dence for Landau levels. Our studies show that for about 
2000 particles, Landau levels will show up as a sequence 
of discernable steps only when 0.99 < Cl/cu± < 1, which is 
achievable with the current capability to control frequen- 
cies, especially for large u>±. For lower frequencies, the 
existence of Landau levels can only be revealed through 
the presence of a "quantized core" , which shrinks in size 
as Cl decreases. On the other hand, the LDA in fig. lb 
shows that by introducing an additional potential V we di, 
Landau levels (in the form of a sequence of steps near 



2 



the center or a "quantized core") can still be revealed 
at frequencies farther beyond lu±, even though the steps 
near the surface are smeared out. 

3D case : For a 3D harmonic trap, eq.(|l|) ac- 
quires terms p 2 /2M + ^Moj z z 2 , which give rise to 
harmonic oscillator eigenfunctions fn z (z) with eigen- 
value e„ 2 = huj z (n z + i). The density is 

P(r,z)= E„,,„, m l/n*(*)| 2 K,m( r )| 2 ®0 - en*,n,m), 

where e„ Z! „. m = e„ z + e n ,m- When £1/uj± ~ 1, we first 
perform the m-sum because it produces the smoothest 
change in the energy. Repeating the steps leading to 
eq.(|TJj), we have 



IfnAzW 



8(/i(r) - hn(2n + 1) - e„J (15) 



where p(r) = p, - Mwx^x - Sl)r 2 w \i- ^{lu 2 ^ - Vt 2 )r 2 . 

The order of summation of the remaining integers 
n z and n depends on the relative strength between to z 
and lo±. When lu z << ujj_, we first sum n z . To do 
that, we note that the density of a ID Fermi gas is 
Pid = \/2Mp/ (irh), where p is the chemical potential. 
In a harmonic trap, the density within LDA is obtained 
by changing p to p(z) = p — ^Mu> 2 z 2 . We then have 

= E„JJU*)I 2 ©(m-0 = ^pl, which 

turns eq.(n3) into a useful LDA 



p{r) 



1 1 



■E 



[/i(r,z) - Ml{2n- 



l)]/(2^). 



p(r) 



71 f l 



|/n,(«)| 2 Int 



/x(r) - huj z (n z + 
2M7 



(17) 



where /ix(r) = /x - \M{u)\ - Q 2 )r 2 - V wa \\. In this limit, 
the density consists of a sequence of discs (labeled by 
rt z ) in the z-direction. Each disc \fn z { z )\ 2 consists of a 
sequence of density steps in the xy-plane reflecting the 
number of filled Landau levels. The behavior of the den- 
sity within each disc in the xy-plane is identical to the 
2D case discussed before. Finally, we note that as tem- 
perature increases, the Landau level structure near the 
surface will first melt away, and the melting will pro- 
ceed toward the center. The temperature below which 
the Landau level effect begin to appear is T — 2Tiu>±/kB, 
which is 3.8 x 1CT 7 K and 9.6 x 1CT 6 K for for to x = 4000Hz 
and 10 5 Hz resp., a temperature range achievable in cur- 
rent experiments ||. 

So far, we have only discussed the effect of Landau lev- 
els on the density profiles of fast rotating Fermi gases. If 
the development of quantum Hall effect in the last decade 
is a guide, one expects many more novel phenomena in 
Fermi gases in the fast rotating regime. This work is 
completed during a workshop at the Lorentz Center of 
University of Leiden. We thank Professor Henk Stoof 
and the Lorentz Center for generous support. This work 
is supported by a Grant from NASA NAG8-1441, and by 
the NSF Grants DMR-9705295 and DMR-9807284. 



(16) 



where a z — y/h/Mu> z , and p(r, z) = p — - £l 2 )r 2 — 

^-uj z z 2 — V wa \\. We have included V wa n in a for the more 
general situation as in the 2D case. Eq.()l6|) shows that 
p(r) is a sum of ID densities (labelled by "n") each of 
which distributed over in a "Landau volume" bounded by 
the "Landau surface" p(r,z) = 7ifi(2n+l). When V^aTi = 
0, SI < oj± the Landau surface are ellipsoidal surfaces. It 
is easy to verify that the surface areas A n for successive 
ellipsoids differ by a constant A„-i—A n = ( ffi-n 2 )' 
When Q > cl>j_, centrifugal instability against harmonic 
confinement sets in and stability can only be established 
by y waU . 

To demonstrate the validity of the LDA eq.jl^), we 
have evaluated the density numerically for a system of 
2000 Fermions for lu z /luj_ = 0.2 at Q/ujj_ = 0.99. The 
results are shown in figure 4. It shows that the LDA 
(dotted line) is a good approximation. The Landau vol- 
umes can be clearly identified by the change of slope in 
the density. The appearance of a plateau the center is 
because u> z /u>± is only 0.2, revealing the 2D feature of 
the n z levels. For smaller ratios of u x /u i , the plateau 
disappears and the LDA expression (eq. ( |16|)) is achieved. 

When u> z > uj±, summation of n in eq.(|l5|) gives 
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The extent of erf(x) is about 1.5, we then have 
iA(s/ v /2mS) = L5 - Since A ( s ) 



2rA(r)/aj_, we have 



\/ m o A ( r ) ~ 75V m o a - 
Note that po(r;L) = Eo~ 
where u nm[ \— (7m!ra!) - 
(7m!m!)- 1/2 9£[(A - w) rr ' 
we have Em=ol Unm ; A l 2 = 



1/2 



Since u n m= u 



tim;A-0 i 



>/2gn 
| 2 /2l 



(with a 



Auj* — | w | 



i). 



£s (vV4)"En=o 



itly, pi(r;mt) = (1 + |V 2 )po(r; 
I v2 + ^ v4 )Po(r;m2), p 3 {r;ml) 

1 V76\ ' 



( v2 a/ 4 )" [po(r- \;L)e x2 ]. Explic- 
it), P2(r;ml) = (1 + 



(1 + fV 2 + ^V 4 + 

38 4 ' S )Po(r;m,3). 
The studies of finite temperature effects will be presented 
elsewhere. At higher temperatures, the density profile is 



that of a Boltzman gas, oc e 
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Caption 

Fig. la and lb: Fig. la and 1(b) corresponds to 17 < u± 
and > lu± resp. The rapid drop at large r is due to the 
strongly repulsive potential V wa \\ . The LDA densities are 
indicated by the steps in thick lines. The integer value of 
(/x(r)+7iQ)/(27iQ), (i.e. I), is related to the index n of the 
intersected Landau level as I = n + 1. And the relation 
(fi(r)+hQ,)/(2hft) = I is equivalent to /x(r) = (2n+l)Sl. 
Figure. 2a: LDA (dotted line) and numerical calculations 
(solid line) of the density for N — 2000 Fermions at 
n/uj ± = 0.996. 

Figure 2b. Density profiles of N = 2000 Fermions at 
0/Vl= 0.993 (crosses), 0.996 (dots), and 0.998 (solid 
line) . 

Figure 3a: LDA (dotted line) and numerical calculations 
(solid line) of the density for TV = 1000 Fermions at 
fl/iu ± = 0.982. 

Figure 3b. Density profiles of N = 2000 Fermions at 
n/uj ± = 0.982 (crosses), 0.989 (dots), and 0.995 (solid 
line) . 

Figure 4. LDA (dotted line) and numerical calculations 
(solid line) of the density for N = 2000 Fermions at 
fi/VL= 0.99 and uj z /uj±= 0.2. 
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N = 2000 




Q. 1 co ± = 0.99 
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